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Abstract 

To describe rich phenomena of dilute dispersions of highly-charged 
colloid particles, it is necessary to modify the long-range electric part of 
the DLVO potential. The screened Coulomb potential of pure repulsion 
obtained from the Helmholtz free energy in the DLVO theory is replaced 
by a new screened electric potential with long-range attractive tail derived 
from the generalized Gibbs free energy which is identified with a total sum 
of chemical potentials. 

Colloidal dispersions display a variety of phenomena [T] such as ordering for- 
mation and phase transitions. Diversity of those phenomena has its origin in a 
subtle characteristic of interaction potential of colloid particles, vifhich changes 
its medium- and long-range properties depending on salts and particle concen- 
trations. The DLVO potential [2l [3] consisting of the van der Waals attraction 
and the screened Coulomb repulsion which is derived from the Helmholz free 
energy F of the system can explain the Schulze-Hardy empirical rule for co- 
agulation of colloid particles. This success of short-range characteristics of the 
colloid interaction had granted it the position of the standard theory of colloid 
dispersions. However, it turns out impossible to describe rich phenomena of 
medium- and long-range nature of colloid dispersions in the DLVO theory. In 
this note, we show that a generalized concept of the Gibbs free energy G is more 
appropriate than the Helmholtz free energy F to describe the chemical as well 
as thermal equilibrium of the colloidal dispersion and prove that a pair potential 
derived from G possesses a favorable property which can explain characteristics 
of the dispersions over all ranges. 

To maintain analytical derivations of thermodynamic functions of colloid 
systems, we accept the mean field theory developed originally by Debye and 
Hiickel [H [5l [6] (DH) to investigate the strong electrolyte. In their theory, the 
electric part of the internal energy E'^^ is calculated from the mean electric 
potential ^{r) which is determined by solving the linearized Poisson-Boltzmann 
(PB) equation. Thanks to homogeneousness of the system, the thermodynamic 
relations 
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can be applied step by step to derive the electric parts of the Helmholtz free 
energy F"^^ and then the Gibbs free energy G^K 

It is also possible to calculate G"^' by summing up the chemical potentials of 
all ions as 



in which Nj is the number of ions of j-th species with valence Zj. Here, we 
introduce the differential operator 



for convenience in later arguments on the colloid interaction. Note that the 
operator D is compatible with the condition of electric charge neutrality 



in the sense that this relation is invariant under its action. 

The DH method has previously been adopted in theories of colloid disper- 
sions. First, Levine and Dube [7] derived a pair potential of particles with a 
medium-range repulsion and long-range attraction from the internal energy E^^ . 
Subsequently, Derjaguin-Landau [2J and Verway-Overbeek [3] utilized Eq. ([T]) to 
derive F^'' from E^'- and obtained a purely repulsive pair potential. Note that 
the Legendre transformation in Eq. ([1]) is applicable because colloid dispersions 
which are in a thermal bath of solvent can be regarded as homogeneous with 
respect to the temperature T. The final step from F'^'^ to G"^' is still very prob- 
lematic for a colloid. The equivalence of two routes to the Gibbs free energy 
by means of the Legendre transformation in Eq. ^ (route L) and the sum of 
the chemical potentials in Eq. 1^ (route C), which holds in the DH theory of 
strong electrolytes, seems not to hold [H [9j [10] for the inhomogeneous colloid 
dispersions. We will carry out below explicit derivations along these two routes 
and prove that, while the route C gives us a favorable pair potential, the route 
L leads us to an unphysical result. 

To examine the electric interaction, let us focus our attention on a dilute 
dispersion of colloid particles with large analytic surface charges. The large 
analytic charges act to attract and to condense counterions near the surfaces 
and the condensed counterions work to screen the analytic charges. Hereafter, 
we call such screened charge the effective charge and the colloid particle with 
the effective charge simply the particle. In the dilute dispersion where particles 
have static configurations, the gasses of small ions are postulated to be in the 
thermal and chemical equilibrium described by the mean electric potential ip{r) 
satisfying the linearized PB equation. 
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For simplicity, we consider a monodisperse dispersion consisting of particles 
with effective valence Z and radius a. The internal energy E'^^ is given by 
integrating the electric energy density over the dispersion as follows: 



(6) 



where Ujir) is the Boltzmann distribution of the small ions and p(r) is the charge 
distribution function of the particle at i?„. Note that the mean electric potential 
is defined except for arbitrary constants. Namely, the theory is invariant under 
the shift "^(r) — ^ ?/;(r) + c with an arbitrary constant c. By applying this 
replacement to Eq. ^ , we obtain the constraint of charge neutrality 

^z,iV,+7VZ = (7) 
i 

where N is the number of particles. Note that, while iV is fixed, Z is assumed 
to be variable in the present postulation. 

Substituting the mean electric potential determined by solving the linearized 
PB equation into Eq. ^ and integrating out the degrees of freedom of the 
electric field and small ions, we find the following expression 

m^n n 

in which Rmn = \Rm — Rn\ is the center to center distance between the m-th 
and n-th particles. The pair potential takes the form [1] 
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where 

z* = z':^^^ (10) 

Ka 

and K is Debye's screening parameter defined by 
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The quantity in Eq. is the self-energy of the n-th particle in the ionic 
atmosphere of the dispersion. 

In parallel with the relation in Eq. ([S]) , it is possible to express the free energy 
F'^' as 

F''=Fo + IY. U^iRrr^n) + E (12) 

in terms of the Helmholtz pair potential and self-energy. As already remarked 
above, the temperature T is a relevant thermodynamic variable for the colloid 
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dispersion. Therefore, we can obtain the Helmhohz free energy F'^' from 
by way of the relation in Eq. ([1]). Thus, the two kinds of pair potentials satisfy 
the differential equation 

U^iR)^U^iR)^T^^Q^ (13) 
which can readily be solved, resulting in 

U^ = ±^le--R (14) 

Here it should be pointed out that the pair potentials in Eqs. © and are 
approximately equivalent to the ones derived, respectively, by Levine-Dube and 
DLVO. 

Now, the problem is to find a generalized concept of the "Gibbs free energy" 
which is suitable to describe the inhomogeneous colloidal dispersions. The free 
energy G"^' can also be expressed as follows: 

= Go + i E U^'iRrnn ) + ^ (15) 
m^n n 

where U^{R) is the Gibbs pair potential to be derived. It must be examined 
which one of the routes L and C lead us to a suitable pair potential. 

To follow route L, the Legendre transformation in Eq. ([2]) has to be applied 
to the free energies F'^'' and G^' . Consequently, the differential equation 

U^iR)^U^iR)-V^^^ (16) 

is derived for the pair potentials U^{R) and U'^{R). To solve this equation, 
we must observe the facts that the equations in Eqs. and are identical 
under the replacement of ( T and U^{R) ) by { V and U^{R) ) and that the 
Debye's screening parameter through which the pair potentials depend on the 
variables T and V is symmetric with respect to these variables. Therefore, the 
Gibbs pair potential satisfying Eq. ([TBI) is proved to be identical with the electric 
pair potential U^{R) given in Eq. ([5]). Namely, the route L of the Legendre 
transformation is destined to result in the unacceptable result C/*^(i?) — U^{R). 
When combined with the van der Waals potential, this pair potential fails to 
explain the Schulze-Hardy rule. 

The source of failure of the route L seems to be traced back to the fact that 
V in Eqs. ^ and dTTI) is not a suitable thermodynamic variable. In Eqs. ([1]), ^ 
and (|lip . the variables T and V appear exactly symmetric manner. However, 
while the colloid dispersion is homogeneous with respect to T, it is not with 
respect to V. For the gas of small ions, the particles form inhomogeneous envi- 
ronment exerting the volume exclusion effect. In the theory of colloid dispersions 
and macroionic systems, the quantity V must be treated, not a simple variable. 
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but an implicit function of configurations of macroionic particles. Careless ar- 
guments along the route L by treating y as a thermodynamic variable must 
be forbidden. It is prudent to renounce Eqs. 1^ and (IT51) including V as the 
variable and to interpret V in Eq. (1111) simply as a given parameter. 

To realize what is another thermodynamic variable apart from T in the 
treatment along route C, let us apply the neutrality conditions in Eqs. ^ and 
([7]) to regions with various sizes in arbitrary positions in system. If the size of 
the region is not very small, the conditions hold approximately well for both 
strong electrolytes and colloid dispersions. So, this argument implies that the 
mean electric charges inside sizable regions can play the role of a thermodynamic 
variable in ionic suspensions. Therefore, we are allowed to generalize the scheme 
along the route C in the DH theory to colloid. 

In case of strong electrolytes, G"^' is given by the sum of the chemical poten- 
tials of ionic species as in Eq. and the operator D is defined in Eq. ^ so as 
to be compatible with the neutrality condition in Eq. ([5]). Therefore, we should 
define the differential operator D for colloid dispersions that is compatible with 
the neutrality condition in Eq. ([T]) by 
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and the "Gibbs free energy" by 
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where the derivative of F'^^ with respect to valence Z should be interpreted as 
the chemical potential of the effective valence Z which works to measure the 
extent of its fluctuation. In this way, the generalized Gibbs free energy for colloid 
dispersions is defined, along route C, as the total sum of the chemical potentials 
of the small ions and the effective valence of the particles. Substitution of the 
expressions in Eqs. (ITil) and (ITSI) into Eq. (IT51) results readily in the pair potential 
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The chemical and thermal processes advance so as to minimize this pair potential 
in dilute dispersions of highly charged colloid particles. Three types of pair 
potentials U^{R), U^{R) and U'^{R) are compared in figure 1. 

The resultant Gibbs pair potential has a medium-range strong repulsion 
and a long-range weak attraction. The repulsive part explains the Schulze- 
Hardy rule [1] for coagulation and the attractive part describes a variety of 
phenomena of ordering formation in dilute dispersions of particles. It is worth- 
while recognizing that the chemical potential of the effective valence plays the 
role of enhancing the magnitude of the medium-range repulsion. As shown in 
figure 2, the pair potential changes its behavior depending sensitively on the 
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concentrations of the particles and added electrolytes through the parameter 
Ka. Evidently, the simple relation 

U^{R) = U^{R) + U^{R) (20) 

holds among the three pair potentials. The "Gibbs pair potential" derived 
along route C shares the good nature of the pair potentials of the internal and 
Helmholtz free energies. 
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Figure 1: Comparison of three types of pair potentials. 
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Figure 2: Behavior of the Gibbs pair potential. 

In this way, improvement of the long-range features of the DLVO theory has 
been achieved without losing the irreplaceable forte in its short-range charac- 
teristics. Guided by observations on the constraint of charge neutrality, we have 
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succeeded in arriving at a generalized concept of the Gibb free energy for colloid 
as the total sum of the chemical potentials of the small ions and the effective 
valence of the particles. It seems reasonable to evaluate this scheme as a mile- 
stone in trials to improve the fundamental concepts of the thermodynamics of 
homogeneous systems so as to be applicable to describe a variety of phenomena 
in intrinsically inhomogeneous systems. Colloid dispersions which provide us 
with definite experimental results will be indispensable and fruitful in studies 
of complex and inhomogeneous ionic systems. 
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